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ABSTRACT 

We present new ultraviolet, optical, and X-ray data on the Phoenix galaxy cluster (SPT-CLJ2344- 
4243). Deep optical imaging reveals previously-undetected filaments of star formation, extending to 
radii of ~50-T00 kpc in multiple directions. Combined UV-optical spectroscopy of the central galaxy 
reveals a massive (2 x 10 9 M Q ), young (~4.5 Myr) population of stars, consistent with a time-averaged 
star formation rate of 610 ± 50 M 0 yr _1 . We report a strong detection of O VI AA1032,1038, which 
appears to originate primarily in shock-heated gas, but may contain a substantial contribution (>1000 
M@ yr -1 ) from the cooling intracluster medium. We confirm the presence of deep X-ray cavities in 
the inner ~10 kpc, which are amongst the most extreme examples of radio-mode feedback detected to 
date, implying jet powers of 2 — 7 x 10 45 erg s -1 . We provide evidence that the AGN inflating these 
cavities may have only recently transitioned from “quasar-mode” to “radio-mode”, and may currently 
be insufficient to completely offset cooling. A model-subtracted residual X-ray image reveals evidence 
for prior episodes of strong radio-mode feedback at radii of ~100 kpc, with extended “ghost” cavities 
indicating a prior epoch of feedback roughly 100 Myr ago. This residual image also exhibits significant 
asymmetry in the inner ~200 kpc (O.I5R500), reminiscent of infalling cool clouds, either due to minor 
mergers or fragmentation of the cooling ICM. Taken together, these data reveal a rapidly evolving 
cool core which is rich with structure (both spatially and in temperature), is subject to a variety of 
highly energetic processes, and yet is cooling rapidly and forming stars along thin, narrow filaments. 
Subject headings: galaxies: active, galaxies: starburst, X-rays: galaxies: clusters, ultraviolet: galaxies 
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1. INTRODUCTION 


The hot intracluster medium (ICM) is the most mas¬ 
sive baryonic component in galaxy clusters, comprising 
~12% of the total cluster mass, or ~10 13 -10 14 M 0 in 
rich clusters. The bulk of the ICM is very low den¬ 
sity (C 10 -3 cm -3 ), and will require ~10 Gyr to cool 
via thermal Bremsstrahlung radiation. In the centers of 
clusters, however, this situation is reversed. The ICM 
in the inner ~100 kpc can reach high enough density 
that the cooling time is short relative to the age of the 
cluster ( t coo i < 1 Gyr). This ought to result in a run¬ 
away cooling flow, fueling a massive, 100-1000 M 0 yr -1 
sta rburst in the central cluster galaxy (for a review, see 
iFabianl 11994 ). However, despite the fact that roughly 
a third of all galaxy clusters have short central cool¬ 
ing times (e.g . . iBauer et all 12005) lYikhlinin et all 120071 
Hudson et al.l 120101 iMcDonaldl 120111 : IMcDonald et al l 
201.3 lil. such massive starbursts are extremely ra re (e.g., 
McNamara et al.1 120061 IMcDonald et all l2012bj h The 
vast majority of clusters which ought to host runaway 


cooling flows show only mild amounts of star formation 
(I Johnstone et all 119871: iMcNamara fe O’Connelll "" 


Allen 19951: Hicks fe Mushotzkv 20051 OT>ea et al.1 


1989; 


2008 


McDonald_e^_aL J2011 L Hoffer et al.ll2012HDonahue et al l 

20151 iMittal et all 1201 5f >. suggesting that much of the 
predicted cooling at high temperature is not, in fact, oc¬ 
curring. 
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By comparing the cooling r ate of intermediate- 
temperature (~10 5 ~ 6 K) gas (e.g., Bregman et al.ll2001b 


Oegerle et al.ll200lt iPeterson et alJlzuua itiregman et ai.l 

20061: Peterson fo Fabianll2006HSanders et al.ll2010L 1201 lb 


20031: Bregman et al.l 


McDonald et al.l I2014af) to the ongoing star formation 


rate in central cluster galaxies, one can quantify what 
fraction of the cooling ICM is able to condense and form 
stars. These studies typically find that only ~1% of the 
predicted cooling flow is u l timatel y converted into stars. 
Recently, iMcDonald et ali (j2014af) showed that this inef¬ 
ficient cooling can be divided into two contributing parts: 
cooling from high (^10' K) to low (~10 K) temperature 
at ~10% efficiency, and an additionally ~10% efficiency 
of converting cold gas into stars. The former term is 
generally referred to as the “cooling flow problem”, and 
suggests that some form of feedback is preventing gas 
from cooling out of the hot phase. 

The most popular solution to the cooling flow problem 
is that “radio-mode” feedback from the central super- 
massive blackhole is offsetti ng radiative lo sses in_the 
ICM (see reviews bv lMcNamara fe Nulsenll2012HFabianl 
I2012D . Radio jets from the active galactic nucleus 
(AGN) in the central cluster galaxy can inflate bubbles 


the hot gas (e.g.. Birzan et al 

r2004[ 2008f Dunn et al. 

20051 Dunn & Fabian 2006 

Raffertv et al. 2006; 

Nulsen et al. 2007: Cavagnolo et al. 20101 Dong et al. 

2010 

: O’Sullivan et al.l 20111: Hlavacek-Larrondo et al. 

2012 

, 2014|). The ubiquity of radio jets in so-called 


two are intimately linked, while the correlation of the 
feedback strength (e.g., radio power, mechanical energy 
in bubbles) with the X-ray cooling luminosity provides 
evidence that this mode of feedbac k is, on averag e, suf¬ 
ficient to fully offset cooling ( e.g., IRaffertv et all[200(1 


IHlavacek-Larrondo et al.|l2012f) . Further, recent studies 
of newly-discovered high-z clusters i n t he So uth Pole 
Telescope 2500 deg 2 survey (iBleem et al.l l2015h suggest 
that this energy balance between ICM cooling and AGN 
feedback has been in place since z ~ 1 (IMcDonald et al.l 
2013b; H lavacek-Larrondo et al.ll2014f) . 

The study of outliers in a population can often provide 
new insights into the physical processes governing said 
population. Amongst cool core clusters, there is perhaps 
no more extre me system than the Phoenix cluster (SPT- 
CLJ2344-4243; IWilliamson et all 1201 It IMcDonald et all 
120121)1) . This massive (M 500 ~ 1.3 x 10 15 M 0 ) sys¬ 
tem is the most X-ray luminous cluster yet discovered 
(L 2 _iokeV,500 = 8.2 x 10 45 erg s _1 ), with a predicted 
cooling rate of ~2000 M 0 yr _1 . However, contrary to the 
norm, the central galaxy in the Phoenix cluster h arbors 
a mas sive starburst (~800 M 0 yr _1 ; IMcDonald et all 
l2013alh a massive reservoir of molecular gas (Mh 2 ~ 
2 x 10 10 M 0 ; IMcDonald et al.ll2014cD. and a dusty type- 2 
q uasar (QSO; McDonald et al.ll2012bt lUeda et all 1 201 31) . 


In lMcDonald et all (I2012bl )7 we argued that the only fea¬ 

sible way to bring such a vast supply of cold gas into the 
center of a ri ch cluster was via a r un away cooling flow. 
However, IHlavacek-Larrondo et al.l (120141) showed that 
radio-mode feedback is operating at a level that ought to 
offset radiative cooling in this system, although whether 
this energy has had time to couple to the ICM is uncer¬ 
tain. The fact that the central AGN is simultaneously 


providing strong radiative (quasar-mode) and mechani¬ 
cal (radio-mode) feedback suggests that it may be in the 
process of transitioning from a QSO to a radio galaxy, 
and that the starburst is being fueled by gas that cooled 
before the first radio outburst. However, this interpreta¬ 
tion hinges on a single, shallow X-ray observation with 
the Chandra X-ray Observatory (10 ks), in which both 
the X-ray nucleus and cavities are detected at low signif¬ 
icance. 

In an effort to provide a more complete picture of this 
system, we have obtained deep far-UV spectroscopy and 
X-ray imaging spectroscopy using the Hubble Space Tele¬ 
scope Cosmic Origins Spectrograph (HST-COS) and the 
Advanced CCD Imaging Spectrometer (ACIS) on the 
Chandra X-ray Observatory , respectively. These data 
provide a detailed picture of the young stellar popula¬ 
tions, the intermediate-temperature gas (Ovi; 10 5 - 5 K), 
and the hot intracluster medium. These new X-ray 
data represent a factor of >10 increase in exposure time 
over previously-published observations. We have also ob¬ 
tained deep optical and radio data via Megacam on the 
Magellan Clay Telescope and the Giant Meterwave Radio 
Telescope (GMRT), respectively. Combined, these data 
will provide a much more detailed view of the complex 
interplay between cooling and feedback in this extreme 
system. We describe the reduction and analysis of these 
new data in §2, along with supporting data at IR and 
optical wavelengths. In §3 we present new results de¬ 
rived from these data, focusing on the central AGN, the 
starburst, and the cluster core. In §4 we discuss the im¬ 
plications of these results and present a coherent picture 
of the physical processes at work in this unique system. 
We conclude in §5 with a summary of this and previous 
work, with a look towards the future. 

Throughout this work, we assume Ho = 70 km s -1 
Mpc -1 , Dm = 0.27, and = 0.73. We assume z = 
0.596 for the Phoenix cluster, which is based o n o ptical 
spectroscopy of th e member galaxies dRuel et al.l (2014 
IBleem et al.ll2015h . 

2. DATA 

Below, we summarize the acquisition and reduction of 
new data used in this study, along with a brief summary 
of supporting data published in previous works. 

2.1. Chandra: X-ray Imaging Spectroscopy 

The Phoenix cluster (SPT-CLJ2344-4243) was ob¬ 
served by Chandra in 2011 (obsid: 13401; PI: Garmire) 
for a total of 11.9 ks as part of a large X-ray survey 
of SPT-selected clusters (PI: Benson). To understand 
in detail this extreme system, we have obtained an ad¬ 
ditional 117.4ks (obsids: 16135, 16545; Pis: McDon¬ 
ald, Garmire), resulting in a combined exposure time 
of 129.3 ks (see Figure [TJ and a total of 88,042 counts 
in the central IMpc and in the energy range 0.5-8.0 
keV. All observations were performed with ACIS-I. The 
data from each observation were individually processed 
in the standard manner, including cleaning and filtering 
for background flares, using CIAO 4.6 and CALDB 4.6.3. 
The X-ray background was modeled using three distinct 
components. First, the re-scaled ACIS blank-sky back¬ 
ground observations provide an adequate representation 
of both the particle-induced and unresolved cosmic X- 
ray background. A second component, modeled as a soft 
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(0.18 keV apec) exces s, accounts for Galactic interstel¬ 
lar medium emission (Markcvitc h et all [20031. Finally, 
a third component, modeled as a hard (absorbed 40keV 
BREMSS) excess, accounts for the fact that a larger-than- 
normal fraction of the CXB may be unresolved in shorter 
exposures. The latter two components were simultane¬ 
ously fit to the on-source spectrum and an off-source 
spectrum extracted from a blank region of the target 
field, at a physical separation of >3Mpc from the clus¬ 
ter center. 



100 kpc 


Fig. 1.— Chandra ACIS-I 0.5—8.0 keV image of the Phoenix clus¬ 
ter, representing a total exposure of 129.3ks. This deep exposure 
reveals a relatively relaxed morphology with the majority of the 
counts being concentrated in the central ~100kpc. 

In order to search for structure in the X-ray surface 
brightness, images were made in several bandpasses in¬ 
cluding a soft (0.5-2.0 keV) and hard (4.0-8.0 keV) 
band, which trace the hot ICM and central AGN, re¬ 
spectively. To search for faint structure in the soft im¬ 
age, we subtract a two-dimensional model comprised of 
three beta functions with a shared center, position an¬ 
gle, and ellipticity. The residual image after this model 
was subtracted was smoothed both adaptively, using 
CSMOOTH 1 , and with fixed-width Gaussians. The for¬ 
mer technique provides the highest-quality image over 
large scales, while the latter has a lower (but non-zero) 
likelihood of artificially creating features out of noise. 

Thermodynamic profiles were measured by extract¬ 
ing on-source spectra in concentric annuli, where we 
strive for ~5000 X-ray counts per annulus. These 
spectra were empirically deprojected using DSDEPROJ 
(|Sanders fe Fabianl 120071 : iRussell et al.l I2008D . Depro¬ 
jected spectra were background subtracted using off- 
source regions (the cluster only occupies one of the four 
ACIS-I chips) and fit over the range 0.5-10.0keV using 
XSPEC f Arnaudi Il996ll . We model the X-ray spectrum 
with a combination of Galactic absorption (phabs) and 
an optically-thin plasma (mekal), freezing the absorbing 

1 http://cxc.harvard.edu/ciao/ahelp/csmooth.html 


column (nn) to the Galactic value dKalberla et al.lf2005h 
and the redshift to z = 0.596. In the inner two annuli, 
we add an additional absorbed powerlaw component to 
account for emission from the central X-ray-bright AGN. 

2.2. HST-COS: Far UV Spectroscopy 

The Cosmic Origins Spectrograph (COS) is a medium- 
to-high resolution UV spectrograph on HST with a 2.5" 
diameter aperture. The broad UV coverage and sensi¬ 
tivity of this instrument make it ideal for simultaneously 
detecting coronal line emission (e.g., Ovi) and the stel¬ 
lar continuum in young stellar populations. We designed 
our observational setup (ID 13456, PI: McDonald) to be 
able to detect O vi AA1032,1038 over the full extent of 
the central starburst. This required two offset pointings, 
as depicted in Figure [2] Following the COS Instrument 
Handbook 2 , we model the throughput across the aper¬ 
ture, yielding a two-dimensional estimate of the through¬ 
put over the extent of the central cluster galaxy. Our 
setup yields >75% throughput over the full extent of the 
starburst, with >50% throughput out to the edge of the 
central galaxy. Figure [2] shows our estimate of the two- 
dimensional throughput map compared to optical and 
near-UV images from HST WFC3-UVIS and COS of the 
central cluster galaxy. 

We use the G160M grating with a central wavelength 
of 1577A, providing coverage from 1386-175lA. This cor¬ 
responds to rest-frame 870-1100A at z = 0.596, which 
covers the Ly/3, O VI AA1032,1038, and HeiiA1085 lines. 

Each of the two pointings were observed for 8 orbits, to¬ 
taling ~20 ks on-source for each spectrum. This exposure 
time was chosen to allow simultaneous modeling of the 
O vi lines and the UV continuum, to sufficient depth to 
rule out a 1% efficient cooling flow in the absence of an 
Ovi detection. 

The observed data were binned to ~0.8A/pix to im¬ 
prove the signal-to-noise. Geocoronal emission and resid¬ 
ual background was subtracted using a deep blank-sky 
spectrum with the same observational setup 3 . The re¬ 
sulting spectra were m odeled using t he latest version 
of Starburst99 (v7.0.0; iLeitherer et all I1999T) 4 . These 
models have sufficient spectral resolution (< lA) over 
1000A < A < 1100A to constrain the age and metallic- 
ity of the starburst in the core of Phoenix. We are re¬ 
stricted to the empirically-calibrated models in this wave¬ 
length range, which have metallicities of (Galactic) 
and Z 0 /7 (LMC/SM C). We use the la t est m odels from 
lEkstrbm et al.l (|2012f) and lGeorgy et al.1 (1201311 , which in¬ 
clude rotation of massive stars. The fits were performed 
by grid searching over a variety of populations (continu¬ 
ous/instantaneous star formation, low/high metallicity), 
with the ages ranging from 0-50 Myr in steps of 0.5 Myr. 

For each grid position, we perform a least-squares fit 
of the model to the data, allowing the normalization, 
reddening, redshift, stellar dispersion, and emission line 
fluxes (Ly/3, O VI AA1032,1038, and HeliA1085) to vary. 

The best-fitting model was chosen to minimize y 2 . The 
best-fitting models are shown in Figures [3] and [4] and 
will be discussed in §3.1 and §4.1. 

2 http://www.stsci.edu/hst/cos/documents/handbooks/current/cos_cover 

3 http://www.stsci.edu/hst/cos/calibration/airglow.html 

4 http://www.stsci.edu/science/starburst99/docs/default.htm 
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Fig. 2.— Left panel: HST WFC3-UVIS image of the central galaxy in the Phoenix cluster, from [McDonald et al.l l[20 1 3ai) . Middle panel: 
HST-COS near-UV image on the same scale, showing that the UV continuum comes from scales larger than the central AGN. Right panel: 
Combined spectroscopic throughput from two pointings of HST-COS. White contours show the HST WFC3-UVIS i-band imaging, which 
is well covered by the COS apertures. 


2.3. MegaCam Optical Imaging 

The Phoenix cluster (SPT-CLJ2344-4243) was ob¬ 
served on 30 August 2013 with the Megacam instrument 
on the Clay-Ma gellan telescope at L as Campanas Ob¬ 
servatory, Chile dMcLeod et al.l ll~998h . Exposures were 
taken in g (4x200 sec), r (4x600 sec + 2x120 sec), and 
i (4x400 sec) passbands in standard operating mode and 
cover the inner 12' (5Mpc) of the cluster. For this work 
we consider only the inner ~0.5 / - the remaining data 
will be published in an upcoming weak lensing study 
(Applegate et al. in prep). Images were processed and 
stacked using the standard Megacam reduction pipeline 
at the Smithsonian Astrophysical Observatory (SAO) 
Telescope Data Center. Color photom etry was ca lib rate d 
using Stellar Locus Regression 1SLR; [High et al.l 1200911 . 
For a more detailed description of th e data reduc tion and 
calibration, the reader is directed to [High et all (12012T) . 

2.4. Giant Meterwave Radio Telescope: 610 MHz 

Giant Metrewave Radio Telescope (GMRT) 610 MHz 
observations of SPT-CL J2344-4243 were taken on June 
14 and 15, 2013. The total on source time resulting from 
these two observing runs was about 10 hrs, with a usable 
bandwidth of 29 MHz. The data were reduced with the 
Astronomical Image Processin g System 5 ( AIPS), Parse l- 
Tongue (iKettenis et, al.ll2006lf and Obit dCottonl 2008), 
following the scheme detailed in llntema et all (20091. 
Reduction steps include flagging of radio frequency in¬ 
terference (RFI), bandpass and gain calibration. Fol¬ 
lowing that, several cycles of self-calibration were car¬ 
ried out to refine the calibration solutions. Direction- 
dependent gain solutions were then obtained towards 
several bright sources within the field of view. The data 
was i maged bv dividing the field of view i nto smaller 
facets (|Perlevlll989HCorriwell fe Perlevlll992f) . This cor¬ 
rects for the non-coplanar nature of the array and the 
direction-dependent calibration. For more deta ils about 
the da ta re duction the reader is refe rred to llntema et all 
(|2009f) and Ivan Weeren et al.l (|2014fh 

2.5. Supporting Data: UV, Optical, IR 

5 http://www.aips.nrao.edu 


We also include in this pa per optical multi-band i mag- 
ing from HST WFC3-UVIS (IMcDonald et al.ll2013alf an d 
optical imaging spectroscopy dMcDonald et al.l 12014x3) . 
For a detailed description of a specific dataset, we di¬ 
rect the reader to the aforementioned papers. Relevant 
features of the various datasets are summarized below. 

Optical imaging was ac quired on the cluster core 
for IMcDonald et all (l2013af) in five WFC3-UVIS bands: 
F225W, F336W, F475W, F625W, and F814W. While 
relatively shallow, these data provided the first resolved 
view of the starburst in the central cluster galaxy. Ad¬ 
ditionally, large-area ground-based data at z-band was 
obtained with MOSAIC-II on the Blanco 4-m teles cope 
(for details, see lSong et al.ll20l2l ; iBleem et al.l 1 2 0151) . 

Optical imaging spectroscopy of the central galaxy was 
obtained f rom Gemini-S GMOS in IFU mode, and pre¬ 
sented in IMcDonald et all ( 20143 b These data span 
rest-frame 3000-6000A, and reveal complex emission-line 
nebulae in and around the central galaxy. The morphol¬ 
ogy of this gas is, for the most part, consistent with the 
near-UV morphology. To allow a direct comparison of 
HST-COS and GMOS spectroscopy, we extract aperture 
spectra from the GMOS IFU data, matching the COS 
throughput as a function of radius within the aperture. 
This allows us to simultaneously model the UV and op¬ 
tical continuum, as well as compare measured emission 
line fluxes in the UV and optical. 

3. RESULTS 

Below, we summarize the main results that emerge 
from the analysis of these deep UV, X-ray, optical, and 
radio data. We defer a discussion of these results in the 
greater context of the cluster’s evolution to §4. 

3.1. Dust and Young Stellar Populations in the BCG 

Using aperture-match ed spectra from both HS T-COS 
and Gemini-S GMOS (IMcDonald et al.l l2014dh com¬ 
bined with broadband photometry spanning the spectral 
gap at 1200-3500A (McDonald et al.f2013af h we can con¬ 
strain the relative contribution of young and old stellar 
populations to the UV-/optical continuum, along with 
the effects of reddening on the spectrum over nearly an 
order of magnitude in wavelength. In Figure [3l we show 
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Fig. 3.— UV-through-optical spectrum of the central ga l axy in the Phoenix cluster. The optical IFU spectrum, from McDonald et al. 
(I2014dh and broadband photometry, from [McDonald et al.l (I2013alh have been aperture matched to the far-UV HST-COS spectrum based 
on Figure 121 Thin c olored curves show the best-fit stellar population model reddened by a variety of Galactic e xtinction mode l s from 
Cardelli et al. (1989), while thick orange an d purple lines s how t he best-fit models using extinction models from ICalzetti et al.l (|T994). 
The featureless, grey extinction model from Calzetti et al. (1994) provides a better match to the spectrum, which appears to lack the 
characteristic 2175A absorption feature. The best-fit stellar population model consists of a highly-reddened, young (4.5 Myr), metal-poor 
starburst. A continuous starburst model provides a qualitatively similar fit, although the strength of the Balmer jump (inset) appears to 
be more consistent with a recently-quenched starburst. 



the full UV+optical spectrum for the combined area of 
the two HST-COS apertures (Figure [2]). We overplot 
a series of best-fittin g model spectra, comparing differ- 
ent reddening models (ICardelli et al.lfl989l : ICalzetti et al.l 
119941 ) and star formation histories. All models include a 
variable-mass old (6 Gyr) stella r pop ula tion . We note 
that the extinction curve from ICardelli et al.l (119891 ) was 
calibrated using observations of stars in our galaxy and 
the Magellanic clouds, while ICalzetti et al.l (119941 ) was 
calibrated on nearby starburst and blue co mpact galax¬ 
ies. We find, as did ICalzetti et alj (I1994H . that a rela¬ 
tively gray, “feature-less” extinction curve provides the 
best fit to the data, suggesting that the size distribution 
of the dust is skewed towards larger grains. The data 
also suggest a lack of the characteristic 2175A bump, 
which is the strongest absorption feature in the interstel¬ 
lar medium but tend s to be missing in starburst galaxies 
(ICalzetti et al.lfl994ll . iFischera fe Dopital (1201 lH suggest 
that high levels of turbulence can both flatten the cur¬ 
vature of the extinction law and wipe out the 2175 A ab¬ 
sorptio n feature. This is consistent withlMcDonal d~et al.l 
(l2014cf) . where we show that the ionized gas has signifi¬ 
cant velocity structure (( a v ) ~ 300 km s^ 1 ) and is con¬ 
sistent with being ionized primarily by sh ocks. _ 

As suming the flatter extinction curve of lCalzetti et al.1 
dnn, we find an overall good fit (Xdof = 1-93) to the 

data spanning rest-frame 800-6000A. The best-fitting 
model (Figure[3l purple curve) is a 4.5Myr-old starburst 
with a total zero-age mass of 2.2 x 10 9 M 0 . This cor¬ 
responds to <1% of the total stellar mass of the cen¬ 
tral galaxy. The fit quality is only slightly reduced 
(Xdof = 1-97) if we assume a continuous star formation 
history, with the best-fit model representing a ~317 M 0 
yr^ 1 over the past >15 Myr (Figure [3] orange curve). 
These two models disagree on the strength of the Balmer 


jump at ~3650A (see inset of Figure [3]), with the data 
preferring the instantaneous star formation model. This 
provides marginal evidence that star formation has been 
quenched in the past Myr. 

When the far-UV spectrum (1000A < \ rest < llOOA) 
is considered on its own, the data prefer a younger pop¬ 
ulation (^1-2 Myr), as we show in Figure 0] It is worth 
noting that the age is only weakly constrained in this 
region of the spectrum, due to the relatively short wave¬ 
length range. There is little difference in stellar popula¬ 
tions and reddening between the northern and southern 
COS apertures. This may be due to the considerable 
overlap between the two regions (see Figure [2J or may 
indicate a uniformity in the starburst. With only two 
options for metallicit y when synthesizing high-resolution 
spect ra at far-UV fekstrom et al.l I2012t IGeorgy et al.1 
120131 see also §2.2), we can only weakly constrain the 
metallicity from this spectrum. We note that the so¬ 
lar metallicity spectrum provides a significantly worse fit 
than the spectrum for the Z = ZqII population, with 
the former having significantly (factor of ~5) stronger 
absorption lines which are not present in the observed 
spectrum. We stress that, despite the extremely high star 
formation rate inferred by the UV spectrum, the young 
population constitutes <1% of the total stellar mass in 
the central galaxy. This is consistent with our picture of 
BCG formation, since this mode of star formation must 
be suppressed to prevent central cluster galaxies from 
growing too massive (and luminous) by z ~ 0. 

In summary, the combined UV-through-optical spec¬ 
trum of the central galaxy in the Phoenix cluster reveals 
a highly-reddened (E(B-V) ~ 0.5), young stellar popula¬ 
tion, consistent with being recently quenched (<5 Myr). 

3.2. Coronal Emission: Ovi XX1032,1038 
















































































6 



1000 1050 1100 

Km [A] 


Fig. 4. — Far-UV spectra for the two HST-COS pointings de¬ 
scribed in §2.2. In each panel, the red curves shows the best-fit 
young stellar population model, including emission from Ly/3, O VI 
AA1032,1038, [Sivi]A1049, and HeiiA1085 (blue) and a blank-sky 
background component (grey). The best-fitting stellar population 
model is described in the upper right of each panel. In general, 
these spectra are consistent with a highly-reddened, young (<5 
Myr) population of stars. 

The emissivity curve of O VI AA1032,1038 peaks at 
~10 5 5 K, making it an excellent probe of gas at in¬ 
termediate temperatures between the “hot” (>10 7 K, 
probed by X-ray emission lines) and “warm” (~10 4 K, 
probed by optical emission lines) phases. Emis¬ 
sion from O v i has been detected in several nearby 
galaxy clusters (lOeeerle et al.ll200lHBregman et al .1120061 
IMcDonald et, al.1 l2014a[ h allowing an independent esti¬ 
mate of the ICM cooling rate. Here, we attempt a simi¬ 
lar analysis to these earlier works on the central ~15kpc 
of the Phoenix cluster (see Figure [2]). 

In Figure [4] we show the red-side UV spectra for our 
two HST-COS pointings. These spectra are modeled us¬ 
ing the best-fitting young stellar population described in 
§2.2 and §3.1, along with emission lines including Ly/3, 
O Vi AA1032,1038, and HeiiA1085. While the model is, 
overall, a good match to the data, there are several fea¬ 
tures present in the data that are not reflected in the 
models. In the case of the emission at ~1050A this 
may be due to time-varying airglow emission. Other 
features may be absent in the models due to the fact 




-2.0 -1.5 -1.0 -0.5 

log 10 (He II 4686 / [0 III] 5007) 


Fig. 5.— High-ionization UV-optical emission line ratios for 
the combined HST-COS apertures shown in Figure J2] Over plot¬ 
ted are model expectations for a pure cooling flow I Ferland et al.l 
1998; Chatzikos ct al. 2015), photoionizatio n from a dust-free A GN 
JGroves et al.1 120041) , and radiative shocks (I Allen et al.l 120081) . In 
red we show the measured line ratios from the combined HST- 
COS and Gemini GMOS spectra shown in Figure [3] The blue and 
green points show the residual when photoionization from young 
stars and a central AGN have been removed, respectively. The 
increasingly-large errorbars represent our uncertainty in the details 
of these photoionization models. In purple we show the expected 
shift if a 5,000 Mq yr —1 cooling flow were subtracted from these 
data, demonstrating the inability of these data to adequately con¬ 
strain the cooling rate of the ICM. 

that they are empirically calibrated and lacking the full 
popu lation of stars include d in synthetic optical spectra 
(see iLeitherer et al.|[2014j) . We note that He n A1085 is 
detected in the northern aperture only , consistent with 
previo usly-published optical IFU data (IMcDonald et al.l 
I2014cf) which shows that the HenA4686 peak is offset 
from the nucleus, due, most likely, to a highly-ionized 
wind. 

We measure a total flux in the O VI AA1032,1038 dou¬ 
blet of fovi = 2.46 d= 0.09 x 10 -14 erg s -1 cm -2 and 
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fovi = 2-56 ± 0.10 x 10” 14 erg s _1 cm -2 for the north¬ 
ern and southern apertures, respectively. The placement 
of these apertures was chosen so that the spectra could 
be added (i.e., the combined throughput never sums to 
>1; see Figure ED, meaning that the total Ovi lumi¬ 
nosity in the central ~15kpc of the Phoenix cluster is 
L ovr =_7.55 ± 0.20 x 10 43 erg s _1 . For comparison, 
iBreeman et al.l (120061 ) found Lovi ~ 6 x 10 40 erg s” 1 
for the Perseus cluster, using UV spectroscopy from the 
FUSE satellite. 

Under the naive assumption that 100% of the 
OviAA1032,1038 emission is a result of gas cooling 
through ~10 5 ' 5 K, we can estimate the ICM cooling rate 
at int ermediate temperatures. Using the latest CLOUDY 
code dFerland et al.lll998t IChatzikos et al.l[2015fl , and as¬ 
suming initial plasma properties matched to the mea¬ 
sured values in the inner ~150 kpc of the cluster (Z = 
0.6 Z Q , kT = 10 keV), we find M 0vi = (L 0v i/1.37 x 10 39 
erg s” 1 ). For the measured luminosity quoted above, this 
corresponds to a cooling rate of 55,000 M 0 yr _1 . How¬ 
ever, there are additional sources of ionization that may 
be dominating the line flux here, specifically photoioniza- 
tion from the central A GN and heating from shocks. In 
iMcDonald et al.l (l2014dl we show that both of these ion¬ 
ization sources are contributing to the line flux in high- 
ionization lines such as [O hi] A5007 and Hen A4863. 

In Figure [5] we show the combined (both apertures) 
aperture-matched (between HST-COS and GMOS) line 
ratios for various sets of high-ionization UV and optical 
emission lines. The observed, extinction-corrected line 
ratios are c onsistent with having origins in either dust- 
f ree AG N (Groves et al.l 120041 ) or fast radiative shocks 
(lAllen et ali 120081 b We subtract the contribution to 
these emission lines from star-forming regions, convert¬ 
ing the UV-derived star formation rates to emission 
li ne flu x es for a range of stellar photoionization models 
l|Kewlev et al.ll200Tlf . The larger uncertainties on the SF- 
corrected line ratios reflects our uncertainty in the ioniza¬ 
tion field of these young stars. These revised line ratios, 
with stellar photoionization contributions removed (Fig¬ 
ure [5J blue points), remain consistent with an AGN or 
shock-heated origin. 

We remove the AGN contribution to the emission lines 
by extracting the optical spectrum from our GMOS 
IFU data in a 1" wide aperture centered on the nu- 
cleus. Ass umin g a r ange of AGN photoionization models 
([Groves et al.ll2004 f to estimate the O VI emission, based 
on the observed He IIA4686 flux, we are able to infer the 
amount of, for example, O VI emission coming from the 
nucleus, without having any actual spatial information 
from the COS spectroscopy. The inferred contribution 
to each line from the AGN was subtracted, allowing us 
to estimate AGN-free line ratios in the optical and UV.. 
The resulting line ratios, which have both stellar and 
AGN photoionization removed (Figure [5l green points), 
rem ain consisten t with radiative shocks with v ~ 300 km 
s _1 (lAllen et al.l I2008T) . The large error bars, resulting 
from our combined uncertainty in the stellar and AGN 
photoionization fields, are c onsistent w ith the full range 
of magnetic fields tested bv lAllen et al.1 (l2008fl . This i s 
consistent with our earlier work (IMcDonald et al.ff2014cD , 
in which we argued that the warm (10 4 K) gas is predom¬ 
inantly heated by radiative shocks, based on multiple op¬ 


tical line ratio diagnostics (e.g., [Oiii]/H/ 3, [On]/[0 hi], 
as well as the gas kinematics. In this earlier work, we 
showed evidence for a high-velocity, highly-ionized plume 
of gas extending north from the central AGN, along the 
same direction as our two COS pointings. Such highly- 
ionized, high-velocity signatures are typically not ob¬ 
served in low-z clusters, where the warm gas tends to 
be only weakly ionized. 

Figure [5] demonstrates that, even after removing ion¬ 
ization contributions from young stars and AGN, the ob¬ 
served high-ionization line ratios are inc onsistent w ith 
a pure cooling flow model f rom cloudy (jFerland et al.1 
119981 IChatzikos et al.1 1201511 . However, we demonstrate 
that including a 5000 M 0 yr _1 cooling flow to the model 
(representing only ~10% of the observed O VI flux, purple 
point in Figure 0 does not significantly change the ob¬ 
served line ratios, which remain consistent with shocks. 
Further, we expec t additio nal ioni zati on from both par¬ 
ticle heating (e.g..[Ferland et al.ll2009lf and mixing (e.g., 
iFabian et al.11201 lh ) The former likely contributes at 
a low level throughout the cluster core. Given that 
the radio flux of the AGN is fairly typical of BCGs 
with significantly less star formation, we expect parti¬ 
cle heating to be negligible compared to other ionization 
sources in the extreme environment of the Phoenix clus¬ 
ter. On the other hand, mixing of hot and cold gas is 
likely contributing significantly given the abundance of 
multiphase gas in this system. Mixing of the hot and 
cold gas ought to result in an intermediate-tem perat ure 
phase , which would likely be near 10 6 K (IFabian et al.1 
120111 b This would lead to substantial Ovi emission for 
relatively low rates of net cooling from the hot to cold 
phase. Thus, we conclude that, as a result of the myriad 
of a dditional high-ionization sources, incl uding shocks 
(e.g.,[McDonald et al.ll2012alh mixin g (e.g.,IFabian et al.1 
1201 lib and particle heating (e.g.- lFerland et al.ll2009lh we 
are unable to constrain the properties of the cooling ICM, 
which is contributing negligibly to the total O VI flux. 
Additional high-ionization lines such as O VII or Fe XIV 
would improve these constraints, perhaps allowing an es¬ 
timate of the amount of gas cooling through ~10 6 K. 

3.3. Star-Forming Filaments on 100kpc Scales 

In IMcDonald et all (l2013af) we presented high-angular- 
resolution broadband imaging of the cluster core, re¬ 
vealing complex filaments of star formation on scales of 
^40 kpc. While providing unmatched angular resolution, 
these data were relatively shallow (~12m), making the 
detection of more extended, diffuse emission at large radii 
challenging. In contrast, the MegaCam data presented, 
while having inferior angular resolution to the HST 
data, are considerably deeper, providing substantially- 
improved sensitivity to faint, extended sources. 

In Figure [6] we show a three-color image which com¬ 
bines two MegaCam bands ( g , r) and a redder band 
(z) from the MOSAIC-II camera on the Blanco 4-m 
telescope. In the cluster core, the new data show 
that t he star-forming filamen ts initially identified with 
HST (IMcDonald et al.1 12013af) extend significantly fur¬ 
ther than previously thought. To the south, a pair of 
filaments extend for ~40 kpc each, while to the north¬ 
west a third filament extends for ~65kpc. The most 
extended filament extends due north from the central 
galaxy for ~100kpc (highlighted in Figure 0. This is 
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Fig. 6. — Optical g,r,z image of the inner region of the Phoenix cluster (SPT-CLJ2344-4243). The central cluster galaxy is located at 
the center of the image. Extending radially from this galaxy are several blue filaments, likely sites of ongoing star formation. We highlight 
the longest filament to the north, which is extended for ~100kpc. In the lower right we s how a zoom-in o f the ce ntral galaxy, based on 
higher-resolution HST data in the same bandpasses. These data were originally presented in I McDonald et al.l (12013aD and reveal a complex, 
filamentary morphology in this starburst galaxy. 


the most extended star-forming filament yet detected in 
a cool core cluster, exceeding the well-studied filaments 
in the nearby Pers eus (60kpc; iConselice et al.l [2 001; 
Canning et al.l [2014lh Abell 179 5 (50kpc; ICowie et al l 
198.1t McDonald fc Veillenxll2009fl and R XJ1532.9+3021 
150 kpc: I'Hlavacek-Larrondo et al l [2013T ) clusters. These 
four filaments, particularly the northern pair, are excep¬ 
tionally straight, similar to the northern filament in the 
Perseus cluster. This morphology has been used to argue 
for relatively lo w tu rbule nce in the core of the Perseus 
clusterfe.g.. iFabian et alil’2008 '). although counter argu¬ 
ments can be made on the basis o f ICM density fluctua¬ 
tions (e.g.. iZhuravleva et al.ll2015h . 

The filamentary emission is observed in both the g and 


r bands for all four filaments, providing preliminary evi¬ 
dence that this is continuum, rather than line, emission. 
At the rest frame of the Phoenix cluster, the g-band 
spans 2400-3300A. This wavelength range is relatively 
free of strong emission lines, with MgnA2798 being the 
only line that may be contributing significantly to the 
flux in this band. This line should be faint in most sce¬ 
narios, with the exception of fast (v > 500 km s -1 ) ra¬ 
diative shocks in a dense (n > 100 cm -3 ) medium. 

The total amount of rest-frame UV emission in the 
outer filaments - thos e that were not included in 
IMcDonald et al.l (|2013af ) - is relatively small. For ex¬ 
ample, the outer 50% of the northern filament (50- 
100 kpc) contributes ^0.5% of the total blue emission, 
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while the thin northwestern filament contributes an addi¬ 
tional ~1%. Thus, the global estimate of the star forma¬ 
tion is not significantly altered by including these deeper 
data. 


3.4. X-ray Surface Brightness Maps 



Fig. 7.— Upper panel: Gaussian smoothed (FWHM = 2") 0.5- 
2.0 keV image of the central ~100kpc of the Phoenix cluster. This 
image clearly shows the pair of cavities ~10kpc to the north and 
south of the X-ray peak, despite the lack of any additional process¬ 
ing. Lower panel: Unsharp masked image of the same region as 
above. This residual image highlights the small-scale structure in 
the inner region of the cluster, showing a pair of highly-significant 
cavities to the north and south and overdense regions to the east 
and west. 


In Figure [7] we show a smoothed 0.5-2.0 keV image 
of the central ~100kpc in the Phoenix cluster. Without 
any additional processing, the p air o f cav ities in the i nner 
~10kpc, reported initially bv iHlavacek-Larrondo et al.l 
(IMj , are evident. An unsharp masked image of the 
core (lower panel of Figure[7]) reveals significant structure 
in the inner ~30 kpc, with overdense regions surrounding 
the pair of cavities. The cavities, located ~10kpc to the 
north and south of the central AGN, are detected at a 
significance of ~25cr in this residual image. The unsharp 
mask technique is only sensitive to substructure on a 
single angular scale, so it is unsurprising that we detect 
no large-scale asymmetries beyond the core region in this 
residual image. 

In order to look for ICM substructure on multiple an¬ 


gular scales, we model the 0.5-2.0 keV surface bright¬ 
ness map with a sum of two beta models and a con¬ 
stant background, using SHERPA 6 . The two beta mod¬ 
els, meant to represent the core and outer ICM, share 
a common center, ellipticity, and position angle. We 
note that, since the X-ray spectrum of the central AGN 
is hig hly obscured ([McDonald et al.l l2012bt lUeda et all 
120131) . there is no central point source detected in this 
soft energy band. During the fit, the cavities (Figure 
0 and point sources were masked to prevent any spuri¬ 
ous positive or negative residuals. The best-fitting model 
was subtracted, with the residual images shown in Fig¬ 
ure [8] This residual image shows a significant amount 
of structure in the central ~150kpc. On scales larger 
than ~200kpc (O.I 5 R 500 ) there is little structure in the 
residual image (see rightmost panel in Figure [8]). If the 
cluster had recently experienced a recent minor merger, 
one would expect to find large-scale spiral-shaped surface 
brightness excesses that would be visible on > 100 kpc 
scales in the residual imag e (e.g.. iRoedieer et all 120111 : 
iPaterno-Mahler eFa l. 2013). For comparison, the large- 
sca le spiral feature in Abell 2029 (IPaterno-Mahler et all 
120131 ) represents an excess surface brightness of 22% at a 
radius of ^150 kpc. At similar radii in the Phoenix clus¬ 
ter, the surface brightness uncertainty in the 0.5-2.0 keV 
image is ~15%, meaning that a similar feature to that 
observed in Abell 2029 would be at the level of the noise. 
Thus, we are unable to rule out a recent minor merger. 
In general, the clu ster appears relax ed - an observation 
corroborated bv lMantz et al.1 ( 2015 ?). who find a marked 
absence of any large-scale “sloshing” features in this sys¬ 
tem based on an independent analysis, labeling it as one 
of the most relaxed clusters in the known Universe. 

Directly to the northwest and southeast of the X-ray 
peak there are a pair of cavities, detected with S/N = 25, 
and consisten t in l ocat ion and size to those reported in 
IHlavacek-Larrondo et alJ (I2014D . The northern cavity 
appears to be “leaking” to the west, while the southern 
cavity is slightly extended to the east. If these cavities 
are indeed not confined to a simple bubble morphology, 
then th e est i mates of the AGN power derived in previous 
works iHlavacek-Larrondo et alJ 120141 ) are significantly 
underestimating the full power output of the AGN. In¬ 
terestingly, the four most extended, linear filaments of 
young stars (Figure [6]) appear to trace the AGN outflow, 
with the northern and southern cavities occupying the 
space between the two filaments in the respective direc¬ 
tion. This may be due to the expansion of the cavity 
compressing the surrounding gas, leading to more rapid 
cooling along the edges of these bubbles. 

We note that these two X-ray cavities are at a com¬ 
mon distance of ~18kpc from the X-ray peak, which is 
relatively small given the massive size of the cool core 
(~100-200 kpc). This may provide the first clue towards 
understanding the extreme level of star formation in the 
core of the Phoenix cluster. Under the assumption that 
these bubbles rise buoyantly, the distance of the cavity 
from the radio source can be linked to the age of the 
AGN outburst. This would signal that the outburst of 
AGN feedback in the core of the Phoenix cluster is rela¬ 
tively recent (^10-100 Myr), and may not yet have had 
time to distribute the energy required to offset cooling 


http://cxc.harvard.edu/sherpa4.4/ 
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Fig. 8. — Residual X-ray surface brightness images (0.5-2.0 keV) after a double beta model (shared center, cllipticity, position angle) 
has been subtracted. Left: Residual image smoothed with a fine kernel (Gaussian FWHM = 2.2"). This image highlights the negative 
residuals (red) near the central AGN (white cross) which are assumed to be X-ray cavities. The dashed radial lines depict the locations 
of the most extended star-forming filaments (Figure [BJ. This panel demonstrates that the most extended star-forming filaments appear to 
avoid the central X-ray cavities. There are also three regions of enhanced (blue) surface brightness, surrounding the X-ray cavities, which 
may be sites of enhanced cooling. Middle: Residual image smoothed with a medium kernel (Gaussian FWHM = 3.5"). This image is 
zoomed out by a factor of 3.2 to highlight large-scale features, such as the spiral-shaped overdensities (blue) and the asymmetric central 
cavities (red). Right: Residual image smoothed with a coarse kernel (Gaussian FWHM = 7"). This image is zoomed out by an additional 
factor of 2, to show the relati ve lack of structure on scales larger than ~200kpc (0.15Rsoo)- We highlight a potential set of “ghost cavities” 
(e.g.. [McNamara et al.1120011 '1 to the north and south, '~100 kpc from the cluster center. These may be the remnants of a prior epoch of 
strong radio-mode feedback, although they are only marginally detected (see Table [2ll. 


to the surrounding ICM. We will return to this line of 
reasoning in §4.3. 

The residual X-ray images (Figure[8]) also show a signif¬ 
icant amount of excess above a smooth component. We 
note that this excess is not an artifact induced by the 
sharp negative residuals (cavities), as these have been 
masked in the 2D fitting process. The left-most panel 
of Figure [5] shows three distinct regions of excess emis¬ 
sion (colored blue). These overdensities may be gas that 
has been “pushed aside” by the inflating bubbles, in¬ 
falling cool clouds, or low-mass groups that are in the 
midst of merging with the main cluster. Surrounding 
the northern cavity, the largest of these overdensities ex¬ 
tends over >100 kpc (see central panel of Figure [5). The 
morphology of this overdense region is reminiscent of an 
infalling cloud of cool gas. Alternatively, the fact that 
both the northern bubble and the northern overdensity 
share the same curvature suggests that this may sim¬ 
ply be the outer rim of a bubble expanding in a dense, 
cool core. To classify these with more certainty, we re¬ 
quire deep enough X-ray data to produce temperature 
and metallicity maps. 

On larger scales, there may be a pair of “ghost” cavities 
to the north and south of the cluster center. These are 
located at distances of ~100 kpc from the X-ray peak and 
are detected most readily in a heavily-smoothed residual 
image (right panel of Figure [8]). We will return to a 
discussion of the significance of these cavities in §4.3. 

3.5. Radio Jets, Ionized Outflows, and X-ray Cavities 

In iMcDonald et all (l2014d) . we present evidence 
from spatially-resolved optical spectroscopy of a highly- 
ionized outflow to the north of the central cluster galaxy. 
This outflow was identified via the presence of a peak in 



Fig. 9.— This figure shows the smoothed residual maps from Fig¬ 
ure [8] with 600 MHz radio contours overlaid in white. The beam 
shape of the 600 MHz data is shown in the upper right corner. 
These contours highlight radio emission centered on the X-ray nu¬ 
cleus and elongated in the north-south direction, in the same direc¬ 
tion as the most significant set of cavities. At large radii, the radio 
emission appears to bend in the direction of the “leaky” cavities. 
The extended radio emission to the southeast appears to coincide 
with the most extended cavity, at a distance of ~ 100 kpc. We 
al so hi ghlig ht the di rection of the ionized wind (inset), discussed 
in IMcDonald et al.l (12013aD , with a yellow arrow. 

high-ionization optical emission lines (e.g., [Oill], Hen) 
roughly 20 kpc north of the galaxy nucleus. In Figure [9] 
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TABLE 1 

Deprojected Thermodynamic Profiles from X-ray Spectra 


Radii [kpc] 

kT [keV] 

P [keV cm 3 ] 
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Note. — The central 1" has been clipped due to the presence of 
an X-ray point source (central AGN). 


(inset), we show that this highly-ionized gas lies in the 
same direction as the northwestern cavity in the X-ray 
emission. This may indicate that strong AGN feedback, 
which is excavating the central cavities, is also heating 
the multiphase gas, presumably via shocks^ This is sim¬ 
ilar to what is observed in Abell 2052 dBlanton et al.l 

i2on . 

Figure [9] also shows the extended morphology of the 
600 MHz radio data from GMRT. These data show emis¬ 
sion centered on the X-ray peak and central galaxy, ex¬ 
tended in the same north-south direction as the cavi¬ 
ties. Interestingly, the extended emission spreads to the 
southeast and northwest, similar to the “leaking” cavi¬ 
ties shown in Figure [8] The most extended emission in 
the southeast appears to trace the morphology of the po¬ 
tential “ghost” cavities, identified in the previous section 
(see right panel of Figure (HI). This low-significance cav¬ 
ity may, instead, be an extension of the primary southern 
cavity, with the radio emission leaking out to the south¬ 
east in a similar manner to the northern cavity. This 
seems to be the case based on the central panel of Figure 
[8j Alternatively, the presence of isolated cavities at large 
radii may signal two distinct bursts of AGN feedback, 
separated by a relatively small amount of time. We will 
discuss these scenarios, among others, in §4. 

The S-shaped morphology of the radio d ata ma y also 
be du e to bulk motions of the ICM. Mcndvara l et al.l 
(|2012f ) showed that gas motions in the core of clusters 
could bend MHD jets. Given that the over-dense gas 
shares a similar morphology, it may be that both the X- 
ray residuals and the radio morphology can be explained 
by sloshing of the cool core, as we would expect if this 
system had undergone a recent minor merger. 

We note that, given the large, asymmetric beam in 
these radio data, it is challenging to say with any cer¬ 
tainty how well the radio and X-ray morphologies are re¬ 
lated. We await higher angular resolution radio observa¬ 
tions to establish the morphological connection between 
the X-ray cavities and the radio emission. 

3.6. Thermodynamics of the ICM 

In Figure [TO] and Table [l] we show the results of our 
X-ray deprojection analysis (see §2.2). The measured 
temperature, entropy, and pressure profiles are roughly 
consistent with the profiles derived based on a mass¬ 
modeling approach in [McDonald et all (l2013bl h In gen¬ 
eral, the profiles look as expected for a strong cool core, 
with the temperature dropping by a factor of ~3 be¬ 
tween the peak temperature (r ~ 500 kpc, 0.4R 5O o) 



10 100 1000 
r [kpc] 


Fig. 10.— Thermodynamic profiles fr om X-ray spectroscopy . 
Spectra are deproj ected using dsdeproj ( [San ders & Fabiani 120071 : 
I Russell et a l. 2008) and fit with a combination of absorption from 
neutral gas (phabs) and plasma emission (mekal). Uncertainty 
regions are shown with diamonds. For comparison, we show ex¬ 
pectations (where available) based on averages of low- and high-z 
clusters in similar mass ranges. The central temperature, pressure, 
and entropy in Phoenix are at the extreme end for low-z clusters, 
and are even more extreme compared to other clusters at z > 0. 5 
which tend to have less massive cool cores (IMcDonald et al.l2013bO . 
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and the central temperature. The temperature pro¬ 
file is_shmlar^Jn_sha£e to_tfre_awrage_£rofhes_jr|vlow- 
z (lYikhlinin et ahl 120061 iPratt et all 120071 : iBaldi et all 
l2007t Leccardi fc Molendil 2008bH and hi gh-z clusters 
(IBaldi et all 120121 iMcDonald et al]|2014b ). which typi¬ 
cally peak at 0.3-0.4R 5O o- The factor-of-three drop in 
the temperature in the inner region is on the high end 
of what is observed in massiv e clusters (lYikhlinin et al.l 
120061: IMcDonald et alJfeOMbl) . 

The central pressure (1.1 keV cm -3 ) is exceptionally 
high the highest measured in any cluster to date - 
consistent with the expectation for a strong cool core 
embedded in one of the most massive clusters known. 
The pressure drops off rapidly with radius, consistent 
with the expectat ion from the_ uni vers al p ressu r e pro¬ 
file (e . g.JArnaud et al.| l201Ct iPlanck Collaboration et ahl 


120131: IMcDonald et al.1 1201 4 bib The pressure profile in 
the Phoenix cluster appears more consistent with z ~ 0 
clusters than other massive clusters at z ~ 0.6 (see Fig¬ 
ure m ■ 

The deprojected central (r < 30kpc) entropy is 19.2 
keV cm 2 , consiste nt with what is found for nearby cool 
core clusters (e.g., iCavagnolo et, al.l [20091 : iHudson et al.1 
l2010f h The entr opy profile shows no evidence of ex¬ 
cess entropy (e.g.. ICavagnolo et al.ll2009ll in the inner re - 
gion, following the “baseline” profile of lVoit et aTI (I2005H 
from the central bin to the outermost bin at ~1 Mpc. 
There is no evidence of flattening in the outer part of 
the entropy prof ile, as is seen in many low-z and high- 
z systems (e.g., iBautz et al.l 120091: iWalker et al.1 120131: 


Reiprich et al.l 20131 lUrban et al. 12014 IMcDonald et al.1 


2014bh. 


The measured metallicity profile is consistent 
with the declining profile foun d in low-z clus¬ 
ters (e.g.. iDe G randi et, aD 12004 IBaldi et all 120071 : 
iLeccardi fc Molendil l2008aj T We note that, given 
the large uncertainties, the metallicity profile is also 
consistent with a constant value as a function of radius. 

The ratio of the cooling time to the freefall time 
( tcooi/tff ) in the cores of nearby clusters correlates with 
the presence of multiphase gas, presumably signaling a 
link b etween the warm and hot phase s (jMcCourt et al.1 
l2012f) . Following iGaspari et al.l (|20 1 2f ) . we estimate the 
ratio of the cooling time to the freefall time in the core 
of Phoenix using the following equations: 


_ 3 (n e + rii) 
co ° l _ 2 n e mA{T , Z) 


(1) 





( 2 ) 


We estimate g(r) assuming that the core is in hydrostatic 
equilibrium, and th e coo li ng fu nction (A(T, Z)) follow¬ 
ing iSutherland fc Dopital (119931! , assuming solar metal¬ 
licity and rij = 0.92n e . The resulting t coo i/tff profile is 
shown in Figure [TT] The minimum value of t coo i/tff ~ 5 
is reached ~20kpc from the central AGN, which corre¬ 
sponds to the radius within which the bulk o f the star 
format ion is contained (see Figure|5]and McDonal d et al.1 
l2013af) . This minimum value is on par with the minimum 
value of t coo i/tff in some n earby co ol core clusters , which 
can vary from ~4 to ~20 dVoit fc Donahudl2014f) . How- 
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Fig. 11.— Ratio of the cooling time ( t coo i ) to the free-fall time 
(tff) as a function of radius for the Phoenix cluster. The shaded 
grey band represents our uncertainty, which is primarily driven by 
uncertainty in the temperature profile. We highlight t coo i/tf f = 10 
with a horizontal dashed line, which appears t o be the approximate 
threshold for condensation in nearby clus ters (McCourt et al l 20 121 ; 
ISharma et al.l [20121 ; IGaspari et alj|2012D . We also highlight the 
maximum extent of four distinct star-forming filaments (Figure 
© with thin vertical lines, and the average extent with a thick 
vertical line. The range of filament radii is shown as a shaded 
blue region. This figure demonstrates that the condensation of the 
cooling ICM may be fueling star formation out to radii of ~100 kpc 
in this extreme system. 



ever, such low valu es o f t r . n ni/tff are typically reached at 
r < 10 kpc (IGaspari et al.ll2012i i. 

Interestingly, t coo i/tff remains below 10 (the approxi- 
mate threshold for ICM condensation in nearby clusters 
dMcCourt et al.ll2012t ISharma et aPl2012t IGaspari et al.l 

I2012H1 out to a radius of ~60kpc, which is nearly the 
maximum extent of the star-forming filaments shown in 
Figure [U This, combined with the presence of large- 
scale star-forming filaments, suggests that condensation 
of the cooling ICM is important in this system out to 
exceptionally large radii. 

We will return to this discussion of the cluster’s ther¬ 
modynamic properties later, to quantitatively compare 
cooling and feedback processes in the cluster core. 

4. DISCUSSION 

We have presented new X-ray, ultraviolet, optical, and 
radio data on the core of the Phoenix cluster (SPT- 
CLJ2344-4243), providing the most complete picture of 
this extreme system to date. Below, we discuss the 
implications of these new data in the context of pre¬ 
vious observa t ions and interpre t ation s ( McDonal d_gt_al.l 
2012H l20l3at llleda et all 120131 IMcDonald et alJl2014cF 
Hlavacek-Larrondo et al. m- 

4.1. A Revised Estimate of the Central Galaxy Star 
Formation Rate 
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In [McDonald et all (|2013aT) . we calculate the star for¬ 
mation rate (SFR) in the core of the Phoenix cluster 
based on rest-frame near-UV imaging from HST, obtain¬ 
ing a value of 800 M 0 yr -1 . In this earlier work, we 
assume that, in the absence of dust, the slope of the 
UV spectr al ene rgy distribution (SED) should be flat 
dKennicuttl i!998f ). This allowed us to compute an ex¬ 
tinction correction based on the UV color. However, the 
amount of UV extinction and the age of the stellar popu¬ 
lation are degenerate with regards to the slope of the UV 
SED. Ideally, we would like an alternate, age-free, esti¬ 
mate of t he reddening, which we o btained via Balmer line 
ratios in iMcDonald et al.1 (l2014cT ). This reddening map 
showed a peak value of E(B — V) = 0.5 at the central 
AGN, dropping to E(B — V) ~ 0.1 in the most extended 
filaments. Applying this reddening map to the rest-frame 
near-UV (~3000A) imaging from HST, assuming an ex¬ 
tinction curve from lCardelli et ahl (|1989T) with Ry = 3.1 
and a r elation between the UV luminosity and SFR from 
iKennicuttl (119981 1. yields a revised SFR of 690 M 0 yr -1 . 

This estima te o f th e star formation rate may be arti¬ 
ficially high. iCalzetti et al.1 (11994) show that starburst 
galaxies tend to lack the rapid rise in the UV extinc¬ 
tion curve that is characteristic of the Milky Way and 
LMC/SMC extinction curves. This is confirmed in Fig¬ 
ure [3l where we show tha t a shallow, featureless extinc¬ 
tion curve (e.g., ICalzetti et al.1 il994 ) provides a much 
better match to the UV-through-optical spectrum of the 
central galaxy in the Phoenix cluster. Using the best-fit 
instantaneous-burst model from Figure El and averaging 
over the lifetime of the burst (4.5 Myr), we infer a re¬ 
duced star formation rate of 490 M 0 yr -1 . 

As we show in Figure [2j the throughput of the com¬ 
bined HST-COS apertures is only ^100% in the inner 
region of the starburst. Thus, an accurate estimate of 
the total SFR based on these data requires an aperture 
corre ction. Based on r est-fra me near-UV HST imaging 
from IMcDonald et al.1 (120134 , we find that only ~26% 
of the UV continuum emission is missed by the com¬ 
bined HST-COS apertures. Assuming that this missing 
flux originates from regions with similar extinction, we 
would infer an aperture-corrected star formation rate of 
660 M 0 yr -1 . More likely, the extinction in the ex¬ 
tended filaments is lower than in the inner few kpc. 
Assuming a much more conservative extinction in the 
aperture correction ( E(B — V ) = 0.15, represe n ting th e 
minimum value measured in IMcDonald et all ( 2014d )) 
yields SFR = 570 M 0 yr -1 . Assuming that the cor¬ 
rect answer lies somewhere between these extrema, we 
infer that the aperture- and extinction-corrected, time- 
averaged star formation rate for the central galaxy in 
the Phoenix cluster is 610 ± 50 M 0 yr^ 1 , where the 
uncertainty here represents only the uncertainty in the 
extinction and aperture corrections. We note that this 
UV-based estima te is consi stent with the recent IR-based 
estimate from lTozzi et al.1 (|2015t ) of 530 ± 50 M 0 yr -1 , 
suggesting that these different methods may be converg¬ 
ing on the right answer. 

The dominant systematic uncertainty in this estimate 
is how, exactly, the star formation has proceeded. For 
example, an instantaneous burst of average age 2 Myr, 
which is preferred by the UV-only spectral fit (Figure [5]) , 
would have a time-averaged star formation rate of 1200 


M 0 yr -1 before aperture correction. On the other hand, 
if the star formation has been proceeding at a constant 
level for the past 15 Myr, the time-averaged rate can be 
as low as ~300 M 0 yr -1 (Figure EJ- This uncertainty 
could be reduced by obtaining deep, rest-frame near-UV 
(~2000A) spectroscopy. At these wavelengths, the stellar 
population synthesis models are considerably more ma¬ 
ture and well-tested, which should lead to overall better 
fits to the data and a tightening of the allowed parameter 
space. 

The Megacam data presented here, while significantly 
deeper, do not provide additional constraints on the 
global star formation rate, since the bulk of the UV emis¬ 
sion is contained in the inner ~30kpc. Instead, these 
data provide improved constraints on the extent and 
large-scale morphology of the star-forming filaments. 


4.2. The Origin of the Star-Forming Filaments 

In the central ~30kpc, the UV emission is highly 
asymmetric (Figure [6]), indicative of a vigorous, turbu¬ 
lent starburst. However, at large radii the young stars 
are oriented alon g t hin, lin ear filaments, akin to sys- 
tems like Perseus (Conselice et al.lfeOOll) and Abell 1795 


( McDonald fc Veilleuxl 120091 ). In " ^McDonald et all 

( 20124 we demonstrate that the fuel for the observed 
star formation must originate within the cluster core. 
A scenario in which cold gas is brought into the core 
via infalling gas-rich galaxies and/or groups is unfeasible 
for this system, given the extreme ICM density (efficient 
ram-pressure stripping) and amount of gas needed to fuel 
such a starburst (i.e., multiple gas-rich compact groups). 
A popular scenario for the origin of star-forming fila¬ 
ments in nearby cool core clusters involves cool gas being 
drawn from the clust er core in the wake of buoyant radio- 
blown bubbles (e.g., iFabian et, al.l 1200(1 iChurazov et al.1 
2013). Figures El and El demonstrate that the most ex¬ 
tended filaments surround, rather than trail, the most 
significant set of cavities, inconsistent with the uplift sce¬ 
nario. Indeed, the fact that the star-forming filaments 
extend beyond these cavities in radius suggests that the 
cool gas was present before the current epoch of feedback. 
However, this does not rule out a scenario in which these 
filaments were uplifted by a previous episode of feedback, 
with the bubbles responsible for this action having risen 
to such radii that they are undetectable. 

It may be that no additional mechanism is neces¬ 
sary to produce the extended, star-forming filaments 
observed in the core of the P h oenix cluster. Recen t 
wor k bv iMcCourt et all (|2012D . iSharma et al.1 (120121) . 
and lGaspari et alT i 20121) has shown that local thermody¬ 
namic instabilities can develop when the ratio of the cool¬ 
ing time to the free-fall time ( t coo i/tff ) is less than 10. In 
Figure HU we show this ratio as a function of radius for 
the Phoenix cluster, assuming hydrostatic equilibrium in 
the computation oft//. The cooling time is shorter than 
10 times the free-fall time over an unprecedented ~60 kpc 
in radius, suggesting that the star-forming filaments may 
be fueled by local thermodynamic instabilities in the hot 
ICM, which rapidly conde nse and then “ra in” down onto 
the central cluster galaxy (jVoit et al.1 12015 ). Further, we 
expect a map of t coo i/tff to be asymmetric, since the 
gravitational potential is roughly spherically symmetric 
(in the core), while the gas density is not (see Figure El- 
For example, the spiral-shaped overdensity at ~50 kpc 
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north of the cluster center will likely have t coo i/tff <10, 
since the azimuthally-averaged value at that radius is 
tcooi/tff ~ 10 and t coo i oc l/n e . Thus, despite the fact 
that t coo i/tff > 10 at r > 60kpc, there are likely over- 
dense regions in the ICM with t coo i/tff < 10 out to the 
full extent of the star-forming filaments at r ~ 100 kpc. 

If the star-forming filaments did indeed condense out 
of the hot ICM, there ought to be evidence of cooling in 
the X-ray spectrum. Using the latest Chan dra da ta, we 
estimate the classical cooling rate following 1 White et all 
lfl997l) . using the equation: 

M(i) = Lx{i) - + 

h(i) + A cj){i) h(i) + A cj)(i) 

(3) 

where the i and i ! indices refer to given annuli, Lx(i) 
is the bolometric X-ray luminosity in a given annulus, 
h(i) = | kT{i)/jim p is the energy per particle-mass of 
the hot gas for a given annulus, and A c/>(i) is the change 
in gravitational potential across shell i. The first term 
on the right corresponds to the cooling rate in the ab¬ 
sence of any additional heating, while the second term on 
the right is the correction factor to account for the fact 
that infalling gas will simultaneously be heated gravi¬ 
tationally. Carrying out this sum in the inner 100 kpc, 
we estimate a classical (luminosity-based) cooling rate 
of 3300 ± 200 M q yr -1 , con sistent with the estimat e 
of 2700 ± 700 M 0 yr -1 from iMcDonald et al.1 (l2013al ). 
Assuming the star formation rate of 613 M 0 yr -1 from 
§4.1, based on the most realistic reddening and stellar 
population model, this implies that ~20% of the pre¬ 
dicted cooling flow is converted into stars. This number 
is very close to the exp ected star formation efficienc y 
from cool gas (10-15%; IMcDonald et al.ll201lL l2014al) , 
leaving open the possibility that cooling at the hot phase 
may be ~100% efficient. However, the deprojected X- 
ray spectra in the inner 50 kpc show no evidence of a 
cooling flow, with a single-temperature model provid¬ 
ing an excellent fit (% 2 = 128.75 for 128 degrees of free¬ 
dom). Li kewise, recent X'.M M-Newton observations pub¬ 
lished by iTozzi et al.l (12015D find relatively weak cooling 
signatures in the spectrum at temperatures of 0.3-3.0 
keV, with measurements by the MOS (620+ 3 4o Mq yr -1 ) 
and PN (210^^5 M 0 yr 1 instruments yielding values 
significantly lower than the luminosity-based estimate. 
This may be because AGN feedback has recently halted 
cooling at high temperatures (see §4.3), with the star¬ 
forming filaments being fueled by previously-cooled gas. 
Indeed, in Figure [3l we show that the optical spectrum 
in the vicinity of the B aimer break is consistent with a 
recently-quenched starburst. Alternatively, a significant 
amount of cooling can be “hidden” in the X-ray spec¬ 
trum, since all of the strong coolin g lines at low tem¬ 
perat ure (Ovu, Ovm, FeXVH XX: [Peterson fe Fabianl 
I2006T) are redsliifted to < 0.6 keV where systematic un¬ 
certainties in the Chandra effective area corr ection are 
high due to contamination (jO’Dell et al.ll2013 b Further, 
the disagreement between the spectroscopic cooling rates 
deri ved by t he pn and MOS detectors on XMM-Newton 
(ITozzi et al.1 120151) suggest that systematic uncertain¬ 
ties in our understanding of the soft X-ray response of 
these detectors (on both Chandra and XMM-Newton) 


will dominate any estimate of the spectroscopic cooling 
rate. Thus, while there is no direct evidence for cooling 
in the low-resolution X-ray spectrum, we can not rule out 
the hypothesis that the star-forming filaments are being 
fueled by local thermodynamic instabilities in the ICM. 

The Chandra X-ray data presented here offer further 
support for a local fuel supply. The short central cool¬ 
ing time ( t coo i < 10 8 yr) falls below the threshold for 
the o nset of cooling insta bilities in the hot atmosphere 
(e.g.. lRaffertv et al.ll2008f ). It also meetings the t rnn i/iff 
criterion for cooling instabilities (e.g. iMcCourt et al. 
2012 ; iSharma et al.ll2012t iGasnari et alJl2012tlVoit et al. 

2015( 1. However, the level of star formation lies below 
that expected from unimpeded cooling. Phoenix har¬ 
bors a powerful central radio-AGN that is apparently 
reducing the cooling rate by ~80%, operating in a self¬ 
regulating feedback loop. The absence of spectral signa¬ 
tures of coolin g through 10 6 K (i.e., Fex vn, O Vii) as in 
other clusters (iPeterson fe Fabianl 1200611 , is inconsistent 
with simple, isobaric cooling models, indicating that the 
gas may be cooling through other channel s (e.g., mix¬ 
ing) or has recently been quenched (e.g., iLi fc Brvanl 
12014 1. Future observations with high-resolution X-ray 
grating/microcalorimeter spectrometers (e.g., Astro-H) 
will help to break this degeneracy by providing firm 
constraints on the rate of radiative cooling from high 
(~10 8 K) through low (~10 6 K) temperatures for a sam¬ 
ple of nearby, cool core clusters. 

4.3. Mechanical AGN Feedback 

To estimate the enthalpy released by the radio jets, 
E = 4PV, we measured the cavity sizes and surrounding: 
press ures (e.g., iChurazov et al.l 120001 : I McNamara et al.l 
2000). The X-ray cavities are seen in the raw image with 
well-defined, elliptical shapes. We measured their loca¬ 
tions and sizes using the residual image after smoothing 
it with a gaussian kernel (FWHM = 3 pixels). Their ra¬ 
dial centroids, projected major and minor axes, and en¬ 
ergetics are presented in Tableland shown in Figure IT21 
We assume the cavities are ellipsoidal volumes with axes 
perpendicular to the plane of the sky and equal to the 
projected major axes (upper limit) and minor axes (lower 
limit). The mean jet power assumes the cavities rose 
buoyantly to their cu rrent l ocatio ns in the pl ane o f the 
sky f ollowing (|Birzan et al.l l2004t iMcNamara fc Nulsenl 
12001 . 

For projected distances of 17 kpc, we find a rise time 
of ~ 1.6 x 10 7 yr, assuming the bubbles rise at the sound 
speed, to a more likely ~ 2 — 6 x 10 7 yr, a ssuming the 
bubbles rise buoyantly dBirzan et al.1 [2004 1. These fig¬ 
ures depend on the local gravitational acceleration mea¬ 
sured from the X-ray pressure profile, assuming hydro¬ 
static equilibrium. The gas pressure and mean gas den¬ 
sity at the locations of the cavities were found to be 
3 x 10~ 9 erg cm -1 and 0.15 cm^ 3 , respectively. For a to¬ 
tal cavity enthalpy of 4 pV = 1.8 —2.7 x 10 6 ° erg, we find a 
mean jet power of 1.7 — 7.2 x 10 45 erg s _1 . This mechan¬ 
ical power places Phoenix among of the most energetic 
AGN outbursts known. The total radio luminosity (10 
MHz - 10 GHz), assuming a constant spectral index of 
a = —1.35 (IMcDonald et al.l[2014df . is 3.6 x 10 43 erg s” 1 , 
or roughly one percent of the jet power, consiste nt with 
other nearby galaxy clusters (|Birzan et al.ll2008r) . 

The total mechanical power output measured in the 
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TABLE 2 

Properties of X-ray Cavities 


Cavity r [ ,; ] 

r [kpc] 


a [kpc] 

b ["] 

b [kpc] 

pV 1 [10 59 erg] 

t buoy [10 7 yr] 

P cav 2 [10 45 erg s- 1 ] 

S/N 3 

Central Cavities 
North 2.6 

17.3 

1.7 

11.3 

1.3 

8.6 

1.9 - 2.5 

2.1 - 5.1 

0.8 - 2.5 

24.7 

South 2.6 

17.3 

2.2 

14.6 

1.3 

8.6 

2.5 - 4.2 

1.8 - 5.8 

0.9 - 4.7 

26.2 

Total (inner) 






4.4 - 6.7 


1.7 - 7.2 


Potential Ghost Cavities 









North 17.0 

114.6 

6.6 

44.8 

4.6 

31.1 

18 - 26 

8.7- 13 

1.8 - 3.8 

7.4 

South 14.4 

97.2 

8.6 

58.4 

2.9 

19.3 

11 - 34 

6.4- 14 

1.0 - 6.7 

8.1 

Total (outer) 






29 - 60 


2.8 - 10.5 



Note. — 

1 : Range reflects uncertainty in the three-dimensional bubble shape. 

2 : Range reflects combined uncertainty in three-dimensional bubble shape and bubble rise time (assuming buoyant rise). 
3 : Based on residual image (Figure [8]). Ratio of cavity depth to noise level at the same radius. 


core of the Phoenix cluster is ~ 2 — 7 x 10 45 erg s _1 . 
This is a factor o f ~4 t i mes l ess than that quoted by 
iHlavacek-Larrondo et al. ( 201 4b This large difference is, 
for the most part, a result of shallower X-ray data in the 



Fig. 12.— Upper panel: Smoothed, residual 0.5-2.0 keV image of 
the inner ~50 kpc of the Phoenix cluster (see Figure [8}. Contours 
highlight structure in the positive (dashed) and negative (solid) 
residuals. The red cross highlights the position of the central AGN, 
while red lines show the size and shape of the ellipses used to 
determine cavity energetics. Lower panel: Similar to upper panel, 
but zoomed out by a factor of ~5 to show the positions of the 
extended “ghost” cavities. These cavities are also visible in Figure 
0 


earlier study, which resulted in a factor of ~2 change in 
pressure and a factor of ~4 change in local gravitational 
acceleration. The new estimate of the mechanical power 
lies among the highest measured in a cluster core the 
core, but is still a factor of ~2-4 times less than the 
energy required to offset cooling (T C ooZ,r<ioofcpc = 9-6 ± 
0.1 x 10 45 erg s _1 ). Given the extreme classical cooling 
rate in the inner 100 kpc (3300 Mg yr -1 ), even a small 
mismatch in energetics over a short time could lead to 
the observed star formation rate of ~600 M 0 yr -1 in the 
central galaxy. 

In Figures [8] and [12] we show the locations of two po¬ 
tential “ghost” cavities, at radii of ~100kpc. These po¬ 
tential cavities are detected with signal-to-noise of ~7 
and ~8 (Table [2]) and are along the same north-south 
direction as the radio jets. Due to the lower detection 
significance, it is possible that only one (or neither) of 
these cavities are real. However, by computing the total 
cavity power assuming they are real, we place a rough 
upper limit on the total amount of power in previous 
outbursts. Following the same procedures as above, we 
estimate a total cavity power of 2.8 — 10.5 x 10 45 erg s -1 , 
or anywhere from a third to 100% of the cooling lumi¬ 
nosity. This range overlaps well with the energy in the 
ongoing outburst at small radii, suggesting that the me¬ 
chanical power has been roughly constant over the past 
two epochs of feedback, spanning ~100 Myr. 

In summary, there is strong evidence for ongoing 
mechanical-mode feedback in the inner ~20 kpc of the 
Phoenix cluster, and weaker, but still convincing, ev¬ 
idence for a prior epoch of feedback roughly 100 Myr 
ago. These episodes of feedback have highly uncertain 
energetics (factor of a few in P ctH ,), but are consistent 
with the energy required to offset cooling losses. This 
scenario is further supported by the fact that the cooling 
time at r > lOkpc is longer than the apparent duty cycle 
of this AGN ('"TOO Myr), suggesting that perfect balance 
in a given feedback epoch is not necessary at larger radii. 
Future studies with significantly deeper X-ray and low- 
frequency radio data will provide a more clear picture of 
both the duty cycle and energetics in these radio-mode 
outbursts. 


4.4. A Transitioning AGN? 
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Phoenix is unusual in that its high radio mechanical 
power is accompanied by a bright quasar with X-ray pho¬ 
ton luminosity 5.6 x 10 45 erg s _1 . This high X-ray flux 
is relatively stable, with no significant fluctuations ob¬ 
served between the initial Chandra and Suzaku obser¬ 
vations in 2011 and 2012, and the most recent Chan¬ 
dra observations in 2014. Radio AGN with compara¬ 
ble radiative and mechanical powers are uncommon, as 
are radiatively-effic i ent AGN at the centers of clusters 
([Russell et alMfiTil : IQ ’Sullivan et aI1l2012t [Walker et al.1 

I2014H . The Eddington luminosity for a 10 9 M 0 black 
hole is 1.2 x 10 47 erg s _1 , which would place Phoenix at a 
few percent of the Eddington rate, close to the expected 
transition between radio-mechanical AGN and quasars 
dChurazov et al.l 120051 : iRussell et al.l 1201 .'ill . Standard 
AGN theory posits that the form of power output de¬ 
pends on the specific accretion rate: objects accreting 
above a few percent of the Eddington rate are dominated 
by radiation and those below are dominated by mechan¬ 
ical outflows. That Phoenix exhibits the characteristics 
of both a quasar and a radio galaxy, suggesting that it 
may be in transition between the two states. 

Based on their discovery of neutral iron emission 
from the AGN using the Suzaku observatory, llleda et all 
(1201311 have argued that the BCG is a Type 2 quasar with 
a high absorption column density. Ass uming a bo lomet- 
ric correction factor of 130 (Marconi et al.M2004f) . they 
argue that the total unabsorbed AGN power would be 
6 x 10 47 erg s _1 . This figure implies a black hole accre¬ 
tion rate exceeding 100 M 0 yr _1 , which would be the 
highest rate known in a brightest cluster galaxy. Such 
a high bolometric correction implies a high infrared lu¬ 
minosity, where most of the radiation from the buried 
AGN should be escaping. The total infrared luminosity 
emerging from the quasar and surrounding star forma - 
tion is only Lir — 3.7 x 10 46 ([McDonald et al.ll2012bH . 
This luminosity is comparable to both the mechanical 
and X-ray energy fluxes, but lies two orders of magnitude 
below the bolometric luminosity quoted by lUeda et al .1 
( 2013 1. We conclude that the bolometric correction to 
the nuclear X-ray luminosity is at most a factor of a 
fe w to ten , consistent with the lower values reported 
bv lYasudevan fe Fabianl ([200711 for high-luminosity AGN. 
Therefore the mechanical power of Phoenix’s radio jets 
is comparable to the radiation emerging from the quasar, 
suggesting that it may currently be undergoing a transi¬ 
tion from “quasar-mode” to “radio-mode”. 

4.5. The Radio Mini-Halo 

The 610 MHz image of the core of the Phoenix clus¬ 
ter previously revealed the presence of a central com¬ 
pact radio source, associated with the BCG, and dif¬ 
fuse 40 0 — 500 kpc emission s urrounding this central 
source (Ivan Weeren et all I2014H . The diffuse emission 
is classified as a radio mini-halo and, at a redsliift of 
0.596 ±0.002, is the most distant known. The integrated 
flux density of 17 ± 5 mJy for the mini-halo results in 
a 1.4 GHz radio power of (10.4 ± 3.5) x 10 24 W Hz -1 , 
scaling with a spectral index of — 1 . 1 . 

Considering the lifetimes of the radio emitting elec¬ 
trons and the large extent of the mini-halos, a form of in- 
situ cosmic ray (CR) production or re-acceleration is re¬ 
quired. An explanation for the presence of CR electrons 


is that turbulence in the cluster core re- accelerates a pop¬ 
ulatio n of relativistic fossil electrons ([Gitti et al.l I2002L 
120041) . These electrons could, for example, have origi¬ 
nated from the central radio source, associated with the 
BCG. It has been suggested that gas sloshing, induced 
by a minor merger event, could generate the required 
turbulence. This is supported by the observed spatial 
correlation between the morphologies of some mini-halos 
and spiral-like sloshing patterns in the X-r ay g as of other 
nearby clusters dMazzott a fc Giacintuccill 2008l ). Numeri¬ 
cal simulations (e.g. JZuHone et al.ll2013| ) provide further 
support for this model. The X-ray residual maps shown 
in Figure [5] do show some spiral structure reminiscent 
of sloshing excited by a minor merger, suggesting that 
Phoenix cluster may have undergone a minor merger re¬ 
cently. The lack of any large-scale asymmetry (rightmost 
panel of Figure [ 8 ]) implies that any merging system must 
be relatively small compared to the main cluster. 

Due to the mass of this system, the necessary turbu¬ 
lent velocities to accelerate CRs can be achieved by the 
relatively-smooth accretion of low-mass systems, rather 
than a single, recent event. Since the Phoenix cluster 
is amongst the most massive clusters known (M 500 = 
1.3 x 10 15 M 0 ; iMcDonald et al.l[2012bll . and turbulent 
velocities scale with halo mass (based on structure for¬ 
mation theory), it is feasible that the Phoenix cluster 
maintains a high enough turbulent velocity in its core 
to continuously re-accelerate a population of relativis¬ 
tic electrons. Alternative models, which can also ex¬ 
plain the presence of radio mini halos, invoke secondary 
electrons that are produced by co llisions between CR 
protons and thermal pr o tons (e. g.. iPfro mmer fo Enfilinl 
12004 iFuiita et, al.l l2007t iKeshet fc Loebl 1201 Cl) . These 
secondary models are only successful in explaining a 
small range of observed properties (e.g., IZuHone et al.1 

l20lil . 

4.6. Is the Phoenix Cluster Unique? 

As one of the most massive cool core clusters known, 
the Phoenix cluster can be considered, in many ways, a 
“scaled-up”, but otherwise normal, system. For example, 
the MOO kpc extent of the star-forming filaments is ex¬ 
treme, being the most extended cool filaments observed 
in a cool core cluster to date. However, we showed in 
IMcDonald et al.l ( 2011 1 that multiphase gas is typically 
observed to a maximum radius of r coo i in nearby clus¬ 
ters, where r coo i is defined at the radius within which 
the cooling time is less than 3 Gyr. For Perseus and 
Abell 1795, this radius corresponds to ~60kpc, while for 
the Phoenix cluster we measure r coo i = 112 kpc. Thus, 
while extreme, the extent of these star-forming filaments 
are scaling as expected with the overall mass of the clus¬ 
ter, and do not challenge our understanding of cool core 
clusters. The same can be said for the extreme cooling 
luminosity, central gas pressure, and radio luminosity. 

The extreme star formation rate in the core of the 
Phoenix cluster (MOO M 0 yr _1 ) may be pointing to 
something unique about this cluster. This starburst 
accounts for ~ 20 % of the classically-predicted cooling 
flow, suggesting that cooling may be procee ding very ef¬ 
f icient ly in this system. For comparison, I O 'Pea et al.1 
([2008D found, for a sample of nearby clusters with star¬ 
forming BCGs, a typical ratio of the star formation 
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rate to the cooling rate of 1.8+J;g% - roughly an or¬ 
der of magnitude lower than what we observe in the 
Phoenix cluster. However, when the most extreme end 
of this population is considered - clusters harboring star- 
burst BCGs such as A bell 1835 (SFR ~ 200 M 0 yr” 1 ; 
iMcNamara et~aD [2f)C)6h . RX J1504.1-0248 (SFR - 140 
M 0 yr -1 ; iQgrean et al.l 120 lTili , and MACS 1931.8-2634 
(SFR ~ 170 M 0 yr -1 ; lEhlert et al.l 1201 ill - this ratio 
jumps to 17 ± 6%. That is, it appears that there is a tail 
of extreme clusters that are cooling rapidly, fueling effi¬ 
cient star formation in the central galaxy. In this context, 
the Phoenix cluster represents the most extreme of (but 
not distinct from) a subset of rapidly-cooling clusters 

The cooling properties of the ICM are also, as one may 
expect, quite different in the Phoeni x cluster t han in a 
typical galaxy cluster. In Figure 3 of lVoit et al.l (j2015T ) . 
the normalized density and cooling time profiles of ~100 
massive galaxy clusters at 0 < 2 < 1.2 are shown. Of 
these systems, only the Phoenix cluster violates the min¬ 
imum floor of t coo i/tff = 10, and it does so by a substan¬ 
tial margin. The violation of this boundary is thought to 
initiate strong AGN feedback via condensation or pre¬ 
cipitation of the cooling ICM. This is also reflected in 
Figure [lQl where we show that the ICM entropy pro¬ 
file in the Phoenix cluster is consistent with the baseline 
entropy profile, which rep resents t he expe ctation in the 
absence of any feedback (jVoit et al.ll2005ll . at all radii. 
When considering all 165 known galaxy groups and clus¬ 
ters with existing Chandra data and temperature s in the 
range 4 keV < Tx < 15 keV, ICavagnolo et all (|2009H 
did not find a single cluster with an entropy profile as 
steep as Phoenix over the range 10 < r < 1000 kpc. 
Indeed, the only other cluste r kno wn with such a steep 
entropy profile is H1821+643 dWalker et al.l[2014lf . which 
also harbors a cluster-centric QSO, presumably fueled by 
a rapidly-cooling core. The consistency between the en¬ 
tropy profile in the Phoenix cluster and the baseline en¬ 
tropy profile suggests that the ICM has been allowed to 
cool in this unique system with relatively little resistance. 

Perhaps providing the link between these two char¬ 
acteristic aspects of the Phoenix cluster - the ex¬ 
treme starburst and rapidly-cooling ICM - is the cen¬ 
tral AGN, which is unusual in its own way. This 
central AGN is one of only a few known cluster- 


centric QSOs dO’Sullivan et al.1 2012| 

Ueda et al.l 20131: 

IKirk et all 20151 Reynolds et al.: 2014 

) and appears to 


be transitioning from quasar-mode to radio-mode. The 
current radio outburst appears to have recently begun, 
with the radio-blown bubbles appearing small in size and 
confined to the innermost regions of the cluster, while 
a previous outburst may have occurred as recently as 
~100Myr ago. This AGN appears to be heavily influenc¬ 
ing the local environment, highly ionizing the cool gas in 
the central ^10 kpc , driving an ionized out flow north of 
the central galaxy (McDonald c t al.ll2014cD . and inflat¬ 
ing bubbles in the dense ICM which may be re-directing 
the extended, star-forming filaments (see Figure 0. Un¬ 
derstanding the extreme nature of this AGN, specifically 
why it appears to have had transitioned from mechanical¬ 
mode to radiative-mode, and now back to mechanical¬ 
mode, will likely provide the key to understanding the 
other oddities in the core of the Phoenix cluster. 


5. SUMMARY 

In this paper we present a detailed, multi-wavelength 
study of the core of the Phoenix galaxy cluster (SPT- 
CLJ2344-4243). The analysis presented here is based 
on a combination of new and archival data at radio 
(GMRT), optical (HST-WFC3, Magellan Megacam), ul¬ 
traviolet (HST-WFC3, HST-COS), and X-ray ( Chandra- 
ACIS-I) wavelengths. The primary results of this latest 
study are summarized as follows: 

• Complex, star-forming filaments are observed in 
the rest-frame ultraviolet to extend up to ~100kpc 
from the central cluster galaxy in multiple direc¬ 
tions. These newly-detected filaments extend a 
factor of ~2 times further than the previously- 
reported maximum extent of star formation. 

• Modeling the combined UV-optical (rest-frame 
900A-6000A) spectrum of the central cluster 
galaxy reveals a massive (2.2 x 10 9 M 0 ), young 
(~4.5 Myr) stellar population. Based on these 
data, we estimate a time-averaged, extinction- and 
aperture-corrected star formation rate of 610 ± 50 
M 0 yr -1 . We note that this estimate can vary by 
a factor of ~2 in either direction by changing our 
assumptions on the star formation history. 

• The best-fitting dust model is significantly “grayer” 
than that measured for the Milky Way and 
SMC/LMC, consistent with observations for 
nearby starburst galaxies and distant quasars. The 
data show no evidence of the 2175A bump, suggest¬ 
ing that there may be a process either destroying or 
preventing the formation of small grains and those 
contributing to this bump (i.e., PAHs). 

• We detect significant O VI AA1032,1038 emission 
(L ovi = 7.55 ± 0.20 x 10 43 erg s _1 ) in the inner 
^15 kpc of the cluster core. This emission is con¬ 
sistent with having origins primarily in the central 
AGN and in shock-heated gas along a northern ion¬ 
ized outflow. We are unable to put constraints on 
what fraction of this emission may originate from 
the cooling ICM - the data are consistent with both 
a complete lack of cooling and a massive (i.e., 5000 
M 0 yr -1 ) cooling flow. 

• We confirm the presence of strong (S/N ~ 25) X- 
ray cavities in the inner 20 kpc of the cluster core. 
The total mechanical energy in these cavities is 
Peat; = 2- 7 x 10 45 erg s _1 , depending on their in¬ 
trinsic shape, making this one of the most powerful 
outbursts of radio-mode feedback known. The in¬ 
ferred jet power from these cavities is slightly less 
than the cooling luminosity (L CO oi ~ 10 46 erg s _1 ) 
in the inner 100 kpc. 

• We find that the bolometric X-ray luminosity of the 
AGN ( Lx.boi = 5.6 x 10 45 erg s” 1 ) corresponds to a 
few percent of the Eddington luminosity, consistent 
with a recent transition from quasars to radio-mode 
AGN. The similarity of the bolometric luminosity 
and the mechanical power further supports this pic¬ 
ture of an AGN transitioning from a radiatively- 
efficient mode to a mechanical mode. 
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• We find evidence for an additional set of X-ray 
cavities at larger radii (~100 kpc), suggesting that 
there may have a been prior episode of radio-mode 
feedback ~100 Myr ago. Assuming that both po¬ 
tential “ghost” cavities are real, this prior episode 
of feedback may have had jet powers of 2.8 - 10.5 x 
10 45 erg s _1 , which is similar to what is measured 
in the inner set of cavities. This implies a relatively 
constant mechanical output of the central AGN be¬ 
tween bursts, with a duty cycle of ~100 Myr. 

• The azimuthally-averaged cooling time of the 
ICM is shorter than the precipitation threshold 
(■ t C ooi/tff = 10) at radii of <50 kpc, suggesting 
that local thermodynamic instabilities in the hot 
ICM may be fueling both the star formation and 
AGN feedback. Dense substructures extending be¬ 
yond 50 kpc likely have exceeded the precipitation 
threshold out to much larger radii (~100-150kpc). 

• We see significant substructure in the inner 
~200kpc (O.I 5 R 500 ) of the cluster. The spiral 
shape of this structure is reminiscent of infalling 
cool clouds. This may be cooling gas that has been 
redirected by the strong mechanical feedback, in¬ 
falling cool group-sized halos, or sloshing of the cool 
core. The presence of a radio mini-halo supports 
the X-ray observations that the inner core is highly 
turbulent. 


• Outside of the cool core (r > O.I 5 A 500 ), the cluster 
appears relaxed and has thermodynamic proper¬ 
ties typical of other clusters at similar mass and 
redshift. 


While illuminating in many ways, these new data 
leave several questions unanswered. It remains unclear 
how such vigorous star formation is sustained in the 
midst of massive, radio-mode outbursts from the cen¬ 
tral AGN. The absolute estimate of the star formation 
rate in the central galaxy hinges on understanding the 
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star formation history, which would benefit from broader- 
band near-UV spectroscopy, while confirmation that star 
formation is being fueled by cooling of the hot ICM 
awaits high resolution X-ray spectroscopy. Confirmation 
of the extended cavities at large radii requires deeper 
X-ray follow-up, alongside high angular resolution low- 
frequency radio imaging. As the most extreme cool core 
cluster known, such follow-up studies of this systems will 
allow a deeper understanding of the complex and ongoing 
war between AGN feedback and cooling in dense cluster 
cores. 
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